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Summary. We discuss the ideal gas like models of a trading market. The effect 
of savings on the distribution have been thoroughly reviewed. The market with 
fixed saving factors leads to a Gamma-like distribution. In a market with quenched 
random saving factors for its agents we show that the steady state income (m) 
distribution P(m) in the model has a power law tail with Pareto index v equal to 
unity. We also discuss the detailed numerical results on this model. We analyze the 
distribution of mutual money difference and also develop a master equation for the 
time development of P(m). Precise solutions are then obtained in some special cases. 

1 Introduction 

The distribution of wealth among individuals in an economy has been an 
important area of research in economics, for more than a hundred years. Pareto 
first quantified the high-end of the income distribution in a society and 
found it to follow a power-law 

P(m) ~ m^ 1+u \ (1) 

where P gives the normalized number of people with income to, and the 
exponent v is called the Pareto index. 

Considerable investigations with real data during the last ten years re- 
vealed that the tail of the income distribution indeed follows the above men- 
tioned behavior and the value of the Pareto index v is generally seen to vary 
between 1 and 3 |2HI3 [HIS]- It is also known that typically less than 10% of 
the population in any country possesses about 40% of the total wealth of that 
country and they follow the above law. The rest of the low income popula- 
tion, in fact the majority (90% or more), follow a different distribution which 
is debated to be either Gibbs E] or log- normal 0] . 

Much work has been done recently on models of markets, where economic 
(trading) activity is analogous to some scattering process [SI HJ El lUl El 
El EH El EH as m the kinetic theory of gases or liquids. 
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We put our attention to models where introducing a saving propensity 
(or factor) JS] for the agents, a wealth distribution similar to that in the 
real economy can be obtained |5J I12| . Savings do play an important role in 
determining the nature of the wealth distribution in an economy and this 
has already been observed in some recent investigations ■ Two variants of 
the model have been of recent interest; namely, where the agents have the 
same fixed saving factor [Hj, and where the agents have a quenched random 
distribution of saving factors • While the former has been understood to a 
certain extent (see e.g, [2011^), and argued to resemble a gamma distribution 
|21|. attempts to analyze the latter model are still incomplete (see however 
|22|'l. Further numerical studies of time correlations in the model seem to 
indicate even more intriguing features of the model. In this paper, we intend 
to analyze the second market model with randomly distributed saving factor, 
using a master equation type approach similar to kinetic models of condensed 
matter. 

We have studied here numerically a gas model of a trading market. We have 
considered the effect of saving propensity of the traders. The saving propensity 
is assumed to have a randomness. Our observations indicate that Gibbs and 
Pareto distributions fall in the same category and can appear naturally in the 
century-old and well-established kinetic theory of gas ^7]: Gibbs distribution 
for no saving and Pareto distribution for agents with quenched random saving 
propensity. Our model study also indicates the appearance of self-organized 
criticality 24 in the simplest model so far, namely in the kinetic theory of 
gas models, when the stability effect of savings ^B] is incorporated. 

2 Ideal-gas like models 

We consider an ideal-gas model of a closed economic system where total money 
M and total number of agents N is fixed. No production or migration occurs 
and the only economic activity is confined to trading. Each agent i, individual 
or corporate, possess money nii(t) at time t. In any trading, a pair of traders 
i and j randomly exchange their money 00 El) such that their total money 
is (locally) conserved and none end up with negative money (m,i(t) > 0, i.e, 
debt not allowed): 



time (t) changes by one unit after each trading. The steady-state (t — > oo) 
distribution of money is Gibbs one: 



nii(t) + rrij(t) 



rm(t + 1) + mj{t + 1); 



(2) 



P{m) = (1/T) exp(— m/T);T = M/N. 



(3) 



Hence, no matter how uniform or justified the initial distribution is, the 
eventual steady state corresponds to Gibbs distribution where most of the 
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people have got very little money. This follows from the conservation of money 
and additivity of entropy: 

P(m 1 )P(m 2 ) = P(mi + m 2 ). (4) 

This steady state result is quite robust and realistic too! In fact, several vari- 
ations of the trading, and of the 'lattice' (on which the agents can be put and 
each agent trade with its 'lattice neighbors' only), whether compact, fractal or 
small- world like leaves the distribution unchanged. Some other variations 
like random sharing of an amount 2m 2 only (not of mi + m 2 ) when mi > m 2 
(trading at the level of lower economic class in the trade), lead to even drastic 
situation: all the money in the market drifts to one agent and the rest become 
truely pauper PHU], 

2.1 Effect of fixed or uniform savings 

In any trading, savings come naturally |18j . A saving propensity factor A is 
therefore introduced in the same model |Hj (see |Zj for model without savings) , 
where each trader at time t saves a fraction A of its money rrii(t) and trades 
randomly with the rest: 

rriiit + 1) = mi(t) + Am; mj(t + 1) = rrij(t) — Am (5) 

where 

Am = (1 - X)[e{mi(t) + m 3 (t)} - m,(f)], (6) 
e being a random fraction, coming from the stochastic nature of the trading. 




Fig. 1. Steady state money distribution (a) P(m) for the fixed A model, and (b) 
Pf(m) for some specific values of A in the distributed A model. All data are for 
N = 200. Inset of (b) shows scaling behavior of Pf(m). 
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The market (non-interacting at A = and 1) becomes 'interacting' for 
any non- vanishing A(< 1): For fixed A (same for all agents), the steady state 
distribution Pf (m) of money is exponentially decaying on both sides with the 
most-probable money per agent shifting away from m = (for A = 0) to M/N 
as A — > 1 (Fig. 1(a)). This self-organizing feature of the market, induced by 
sheer self-interest of saving by each agent without any global perspective, is 
quite significant as the fraction of paupers decrease with saving fraction A and 
most people end up with some fraction of the average money in the market 
(for A — ► 1, the socialists' dream is achieved with just people's self-interest 
of saving!). Interestingly, self-organisation also occurs in such market models 
when there is restriction in the commodity market Although this fixed 
saving propensity does not give yet the Pareto-like power-law distribution, 
the Markovian nature of the scattering or trading processes (eqn. ijljl) is lost 
and the system becomes co-operative. Indirectly through A, the agents get 
to know (start interacting with) each other and the system co-operatively 
self-organises towards a most-probable distribution {ra p ^ 0). 

This has been understood to a certain extent (see e.g, [2312]), an d argued 
to resemble a gamma distribution and partly explained analytically |22j . 

2.2 Effect of distributed savings 

In a real society or economy, A is a very inhomogeneous parameter: the interest 
of saving varies from person to person. We move a step closer to the real 
situation where saving factor A is widely distributed within the population 

P21G31E1]. 
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m 

Fig. 2. Steady state money distribution P(m) for the distributed A model with 
< A < 1 for a system of TV = 1000 agents. The x~ 2 is a guide to the observed 
power-law, with 1 + v — 2. 
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The evolution of money in such a trading can be written as: 

m l (t + 1) = X t m l {t) + ey [(1 - \i)m t {t) + (1 - Xj)mj(t)} , (7) 

mj(t + 1) = \ 3 m 3 {t) + (1 - e y ) [(1 - X^m^t) + (1 - A^m^t)] (8) 
One again follows the same trading rules as before, except that 

Am = ey(l - \j)rrij(t) - (1 - A,)(l - ey )«»<(«) (9) 

here; Ai and Xj being the saving propensities of agents i and j. The agents have 
fixed (over time) saving propensities, distributed independently, randomly and 
uniformly (white) within an interval to 1 agent i saves a random fraction 
Ai (0 < A, < 1) and this Ai value is quenched for each agent (Ai are indepen- 
dent of trading or t). Starting with an arbitrary initial (uniform or random) 
distribution of money among the agents, the market evolves with the trad- 
ings. At each time, two agents are randomly selected and the money exchange 
among them occurs, following the above mentioned scheme. We check for the 
steady state, by looking at the stability of the money distribution in succes- 
sive Monte Carlo steps t (we define one Monte Carlo time step as N pairwise 
interations). Eventually, after a typical relaxation time (~ 10 6 for N = 1000 
and uniformly distributed A) dependent on N and the distribution of A, the 
money distribution becomes stationary. After this, we average the money dis- 
tribution over ~ 10 3 time steps. Finally we take configurational average over 
~ 10 5 realizations of the A distribution to get the money distribution P{m). 
It is found to follow a strict power-law decay. This decay fits to Pareto law 
(JTJ with v = 1.01 ± 0.02 (Fig. 2). Note, for finite size N of the market, the 
distribution has a narrow initial growth upto a most-probable value m p after 
which it falls off with a power-law tail for several decades. This Pareto law 
(with v ~ 1) covers the entire range in m of the distribution P(m) in the limit 
N — > co. We checked that this power law is extremely robust: apart from the 
uniform A distribution used in the simulations in Fig. 2, we also checked the 
results for a distribution 

p(X) ~ |A - A| a , A 07 M, 0<A<1, (10) 

of quenched A values among the agents. The Pareto law with v = 1 is universal 
for all a. The data in Fig. 2 corresponds to Ao = 0, a = 0. For negative a 
values, however, we get an initial (small m) Gibbs-like decay in P(m) (see 
Fig. 3). 

In case of uniform distribution of saving propensity A (0 < A < 1), the 
individual money distribution Pf(m) for agents with any particular A value, 
although differs considerably, remains non-monotonic: similar to that for fixed 
A market with m p (A) shifting with A (see Fig. 1). Few subtle points may be 
noted though: while for fixed A the m p (X) were all less than of the order of 
unity (Fig. 1(a)), for distributed A case m p (X) can be considerably larger and 
can approach to the order of N for large A (see Fig. 1(b)). The other important 
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Fig. 3. Steady state money distribution P(m) in the model with for a system of 
N = 100 agents with A distributed as p(\) ~ \ a , with different values of a. In all 
cases, agents play with average money per agent M/N — 1. 

difference is in the scaling behavior of P/(m), as shown in the inset of Fig. 
1(b). In the distributed A ensemble, Pf(m) appears to have a very simple 
scaling: 

PjW-fl-AJfHl-A)), (11) 

for A — ► 1, where the scaling function T(x) has non-monotonic variation 
in x. The fixed (same for all agents) A income distribution Pf(m) do not 
have any such comparative scaling property. It may be noted that a small 
difference exists between the ensembles considered in Fig 1(a) and 1(b): while 
/ mPf(m)dm = M (independent of A), J mPf(m)dm is not a constant and 
infact approaches to order of M as A — > 1. There is also a marked qualitative 
difference in fluctuations (see Fig. 4): while for fixed A, the fluctuations in time 
(around the most-probable value) in the individuals' money rrii(t) gradually 
decreases with increasing A, for quenched distribution of A, the trend gets 
reversed (see Fig. 4). 

We now investigate on the range of distribution of the saving propensities 
in a certain interval a < A, < b, where, < a < b < 1. For uniform distribution 
within the range, we observe the appearance of the same power law in the 
distribution but for a narrower region. As may be seen from Fig. 5, as a — > b, 
the power-law behavior is seen for values a or b approaching more and more 
towards unity: For the same width of the interval \b — a\, one gets power-law 
(with same v) when b — ► 1. This indicates, for fixed A, A = corresponds to 
Gibbs distribution, and one gets Pareto law when A has got non-zero width of 
its distribution extending upto A = 1. This of course indicates a crucial role 
of these high saving propensity agents: the power law behavior is truely valid 
upto the asymptotic limit if A = 1 is included. Indeed, had we assumed A = 1 
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Fig. 4. 

(b), (c) 
(e), (f). 



1000C 



Time variation of the money of the ith trader: For fixed A market - (a), 
and for agents with specific values of A in the distributed A market - (d), 



in I|1Q(I . the Pareto exponent v immediately switches over to v = 1 + a. Of 
course, Ao f= 1 in (|l(Jfl leads to the universality of the Pareto distribution with 
v = 1 (irrespective of Ao and a). Indeed this can be easily rationalised from 
the scaling behavior (fTTJl : P(m) ~ f Q Pf(m)p(X)dX ~ m~ 2 for p(X) given by 
(|1U|> and m _ ( 2+Q ) if Ao = 1 in (|10|l (for large m values). 

These model income distributions P(m) compare very well with the wealth 
distributions of various countries: Data suggests Gibbs like distribution in the 
low-income range (more than 90% of the population) and Pareto-like in the 
high-income range PJ (less than 10% of the population) of various countries. In 
fact, we compared one model simulation of the market with saving propensity 
of the agents distributed following (fTU|) . with Ao = and a = —0.7 12 . The 
qualitative resemblance of the model income distribution with the real data 
for Japan and USA in recent years is quite intriguing. In fact, for negative 
a values in Ijl0|l . the density of traders with low saving propensity is higher 
and since A = ensemble yields Gibbs-like income distribution ©, we see 
an initial Gibbs-like distribution which crosses over to Pareto distribution 
with v = 1.0 for large m values. The position of the crossover point depends 
on the value of a. It is important to note that any distribution of A near 
A = 1, of finite width, eventually gives Pareto law for large m limit. The same 
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Fig. 5. Steady state money distribution in cases when the saving propensity A is 
distributed uniformly within a range of values: (a) width of A distribution is 0.5, 
money distribution shows power law for 0.5 < A < 1.0; (a) width of A distribution is 
0.2, money distribution shows power law for 0.7 < A < 0.9. The power law exponent 
is v ~ 1 in all cases. All data shown here are for N = 100, M/N = 1. 

kind of crossover behavior (from Gibbs to Pareto) can also be reproduced in a 
model market of mixed agents where A = for a finite fraction of population 
and A is distributed uniformly over a finite range near A = 1 for the rest of 
the population. 

We even considered annealed randomness in the saving propensity A: here 
Aj for any agent i changes from one value to another within the range < 
Aj < 1, after each trading. Numerical studies for this annealed model did not 
show any power law behavior for P(m); rather it again becomes exponentially 
decaying on both sides of a most-probable value. 

3 Dynamics of money exchange 

We will now investigate the steady state distribution of money resulting from 
the above two equations representing the trading and money dynamics. We 
will now solve the dynamics of money distribution in two limits. In one case, 
we study the evolution of the mutual money difference among the agents and 
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look for a self-consistent equation for its steady state distribution. In the other 
case, we develop a master equation for the money distribution function. 




m 

Fig. 6. Steady state money distribution P(m) against m in a numerical simulation 
of a market with N = 200, following equations J7J and @ with = 1/2. The 
dotted lines correspond to m~( 1+v ^; v = 1. 



3.1 Distribution of money difference 

Clearly in the process as considered above, the total money (m, + rnj) of the 
pair of agents i and j remains constant, while the difference Arriij evolves as 

(Amij)t+i = ( m i ~ ™>j)t+i 

= (^^) (^)* + (^T^) ^ + m ^ 

+(2ey - 1)[(1 - \i) mi (t) + (1 - X^rrijit)}. (12) 

Numerically, as shown in Fig. 2, we observe that the steady state money 
distribution in the market becomes a power law, following such tradings when 
the saving factor Ai of the agents remain constant over time but varies from 
agent to agent widely. As shown in the numerical simulation results for P(m) 
in Fig. 6, the law, as well as the exponent, remains unchanged even when 
€ij = 1/2 for every trading. This can be justified by the earlier numerical 
observation [H1E1 for fixed A market (A-j = A for all i) that in the steady state, 
criticality occurs as A — > 1 where of course the dynamics becomes extremely 
slow. In other words, after the steady state is realized, the third term in H2|) 
becomes unimportant for the critical behavior. We therefore concentrate on 
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this case, where the above evolution equation for Am^ can be written in a 
more simplified form as 

(Am,ij)t+i = atij(Amij)t + Pij(rm + mj) t , (13) 

where = |(A, + Xj) and fiij — — Xj). As such, < a < 1 and 

-|</3<|. 

The steady state probability distribution D for the modulus A — \Am\ of 
the mutual money difference between any two agents in the market can be 
obtained from 1(13(1 in the following way provided A is very much larger than 
the average money per agent = M/N. This is because, large A can appear 
from 'scattering' involving mi — rrij = ±A and when either rrii or rrij is small. 
When both rrii and rrij are large, maintaining a large A between them, their 
probability is much smaller and hence their contribution. Then if, say, m, is 
large and rrij is not, the right hand side of 1)13(1 becomes ~ (o^ + f3ij)(Aij) t 
and so on. Consequently for large A the distribution D satisfies 

D(A) = [ dA' D(A') (6(A - (a + (3) A 1 ) + 5(A -(a- /3)A')) 




where we have used the symmetry of the (3 distribution and the relation 
otij + ftij = Ai, and have suppressed labels i, j. Here (...) denote average over 
A distribution in the market. Taking now a uniform random distribution of 
the saving factor A, p(A) = 1 for < A < 1, and assuming D(A) ~ z\~( 1+ t) 
for large A, we get 

1 = 2 J dAA 7 = 2(l + 7 )- 1 , (15) 

giving 7 = 1. No other value fits the above equation. This also indicates that 
the money distribution P(m) in the market also follows a similar power law 
variation, P(m) ~ m~( 1+ ^ and v = 7. We will now show in a more rigorous 
way that indeed the only stable solution corresponds to v = 1, as observed 
numerically pilT31IH|. 



3.2 Master equation and its analysis 

We also develop a Boltzmann-like master equation for the time development 
of P(m,t), the probability distribution of money in the market [251 126] . We 
again consider the case = ^ in J7J and JSJ) and rewrite them as 

Collecting the contributions from terms scattering in and subtracting those 
scattering out, we can write the master equation for P(m, t) as 
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dP(m,t) f f 

— \-P(m,t) — ( drrii I dnij P(mi,t)P(mj,t) S(^fmi+fxjmj—m)), 

(17) 

which in the steady state gives 



P{m) = ( J dnii J drrij P(rrii)P{mj) 8([ifmi + fij mj — to)}. (18) 
Assuming, P{m) ~ m~( 1+ly ) for m — > oo, we get [251 I2(i | 

1 = {{ii+y + oo") = 11 dv + d^p(v + Wp-) [(m + ) v + (a* - )"] • (is) 



Considering now the dominant terms (oc x _r for r > 0, or oc ln(l/a;) for r = 0) 
in the x — > limit of the integral J TO^ +r 'P(m) exp(— mx)dm, we get from 
eqn. Ijl9|l . after integrations, 1 = 2/(^+1), giving finally v = 1 (details in 
Appendix). 



4 Summary and Discussions 

We have numerically simulated here ideal-gas like models of trading markets, 
where each agent is identified with a gas molecule and each trading as an 
elastic or money-conserving two-body collision. Unlike in the ideal gas, we 
introduce (quenched) saving propensity of the agents, distributed widely be- 
tween the agents (0 < A < 1). For quenched random variation of A among the 
agents the system remarkably self-organizes to a critical Pareto distribution 

of money with v ~ 1.0 (Fig. 2). The exponent is quite robust: for savings 
distribution p(X) ~ | A — A| a , Ao ^ 1, one gets the same Pareto law with v = 1 
(independent of Ao or a). 

A master equation for P(m,t), as in (|17JI . for the original case (eqns. Q 
and JHJ) was first formulated for fixed A (Aj same for all i), in [23j and solved 
numerically. Later, a generalized master equation for the same, where A is 
distributed, was formulated and solved in [22] and We show here that 
our analytic study clearly support the power-law for P(m) with the exponent 
value v = 1 universally, as observed numerically earlier [1211131 [Li] . 

It may be noted that the trading market model we have talked about here 
has got some apparent limitations. The stochastic nature of trading assumed 
here in the trading market, through the random fraction e in ©, is of course 
not very straightforward as agents apparently go for trading with some def- 
inite purpose (utility maximization of both money and commodity). We are 
however, looking only at the money transactions between the traders. In this 
sense, the income distribution we study here essentially corresponds to 'paper 
money', and not the 'real wealth'. However, even taking money and commod- 
ity together, one can argue (see ^U]) f° r the same stochastic nature of the 
tradings, due to the absence of 'just pricing' and the effects of bargains in the 
market. 
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Apart from the observation that Gibbs (1901) and Pareto (1897) distri- 
butions fall in the same category and can appear naturally in the century-old 
and well-established kinetic theory of gas, that this model study indicates 
the appearance of self-organized criticality in the simplest (gas) model so far, 
when the stability effect of savings incorporated, is remarkable. 
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A Alternative solution of the steady state master 
equation (|18|) 



Let S r (x) = dmP(m)m u+r exp(-mx); r > 0, x > 0. If P(m) = A/ 



m 



l+u 



then 



/>OC 

S r (x) = A / dm m r ~ x exp(— mi) 
Jo 



A 



if r > 



Aln ( - 1 if r = 0. 

x 



(20) 



From eqn. I|18(l . we can write 

S r (x) — 



o 

oo 



(/ drrii / drrij P(mi)P(mj)(rriifif + TOj/i- ) r exp[— (771^/1^ + m^/x - )x}) 



drrii Am\ 1 (exp(-TO i /j+a;) (/j+)' y+r ) 



+ / drrij Am r j 1 (exY>{-mjiL J x) (/x - ) v+r ) 



or, 



drrij P(rnj){exjp(—mj^j x)) 



drrii P(rrii){exp(-mi[if x)) 



(21) 



S r (x) = / di4p{l4)\ I drrii Ami 1 exp(-mi/x+^) 



\J0 



+ / dii j g(/ij ) ( / drrij Arrfj 1 exp(-m j /x j x) ) (/^ ) ""^ (22) 
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since for small x, the terms in the square brackets in 1|21|) approach unity. We 
can therefore rewrite (1221 as 



S r {x) = 2 



dn + (fi + y +r S r (xn+) + / dpT {n-) v+r S r {xn~) 



(23) 



Using now the forms of S r (x) as in 1)20(1 . and collecting terms of order x r 
(for r > 0) or of order ln(l/a;) (for r = 0) from both sides of l(23|l . we get l(19|l . 
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